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ABSTRACT 



This study demonstrates the significance of water vapors 
influence on satellite-retrieved aerosol characteristics using 
NOAA’s AVHRR. An improvement to optical depth and Aerosol 
Particle Size Index (Si 2 ) estimations derived from channels 1 
(0.63pm) and 2 (0.86pm), is made through knowledge of column 
water vapor derived from channels 4 (10. 8pm) and 5 (12.0pm). 
A 2.0 gm cm'2 column water vapor produces a 5% increase in 
S 12 . This results in a 15% error in the variable scattering phase 
function P(0) and retrieved aerosol optical depth (5 a ). The error 
introduced by water vapor is quantified through use of the 
LOWTRAN7 atmospheric propagation model to be applied as a 
reformulated parameterization of (P(@)). 
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I. INTRODUCTION 



A. PURPOSE 

Direct measurements of aerosol characteristics can be provided by ship or 
aircraft. Although they provide valuable information, they are very expensive, 
require detailed coordination, and only small sections of the world's oceans can 
be monitored at any time. Satellites can provide near continuous spatial and 
temporal data desirable for evaluating large scale global effects as well as timely 
local measurements suitable for military applications. Meteorological satellites, 
however, have been designed to measure cloud albedos while much lower 
radiances are associated with aerosol optical depth variations (Durkee 1984). 
Therefore, careful processing of satellite data is necessary to determine aerosol 
characteristics. 

Frost (1988) has developed an algorithm to derive aerosol particle 
characteristics from the Advanced Very High Resolution Radiometer (AVHRR) 
on the NO A A polar orbiting satellites. By taking advantage of the varying 
sensitivities of aerosol scattering in two spectral windows, a multi-channel 
comparison provides information about the aerosol particle size distribution. By 
applying this information to a variable scattering phase function, improved 
estimations of optical depth are possible. 

This study will provide an improvement to Frost’s technique. Frost requires 
analyzed pixels to be free of clouds. However, these clear pixels still contain 
varying amounts of water vapor. Frost assumes that varying amounts of water 
vapor do not appreciably effect retrieved radiances. This study will use the U. S. 
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Air Force’s LOWTRAN7 atmospheric propagation model to demonstrate the 
degree that water vapor influences present aerosol estimations. 

This thesis will: 

1 . Show the effect of varying water vapor amounts on satellite- 
retrieved aerosol characteristics, specifically the particle size index 
parameter or S 12 as defined by Frost (1988). 

2. Provide a parameterization of S 12 corrected for column 
water vapor (w), using a multichannel analysis of AVHRR data. 

B. MOTIVATION 

Radiation in the atmosphere is attenuated by many processes. It is absorbed 
and scattered by molecules, aerosols, and meteorological features such as clouds, 
fog, and precipitation. Some of these influences may be of short duration and of 
local scale such as a cumulus cloud passing overhead and casting a shadow. 
Others, such as stratospheric dust particles are semipermanent and of global 
proportions. 

Aerosols are solid or liquid particles dispersed in the atmosphere. Dust, 
smoke, and even fog or rain are examples of aerosols. Aerosol particles are 
present in the atmosphere in various forms and concentrations. In this thesis, the 
change in particle size due to coalescence of water vapor will not be considered. 

With the ever increasing use of electro-optical devices, military applications 
benefit significantly from improved aerosol measurements. The optimum 
waveband for a laser is from 0.4 to 1.0 pm. This waveband avoids ozone 
absorption at shorter wavelengths, as well as the water vapor absorption that 
dominates longer wavelengths (Bloembergen, et al 1987). This waveband, 
however, is significantly influenced by aerosol scattering. 
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Another important hypothesis that demonstrates the need for better aerosol 
estimations and modeling on a much larger scale is that of Charlson, et al. 
(1987). Charlson proposes that the major source of cloud condensation nuclei 
(CCN) over remote unpolluted oceans appears to be aqueous dimethylsulphide 
gas (DMS) produced by phytoplankton as a waste product. DMS is emitted 
through the sea surface through a chain of complex biological and chemical 
processes where it is oxidized in the atmosphere to form sulfate aerosol particles 
which act as CCN. Increased phytoplankton populations and the resulting 
increase in aerosol particles effectively cool the atmosphere by increasing 
backscatter to space. A 30% increase in aerosols over the ocean surfaces 
worldwide can conservatively decrease the average global temperature by 1.3 K 
(Charlson 1987). 

For many applications, aerosols result in noise interfering with 
measurements of other processes. These applications benefit from improved 
aerosol estimations that eliminate aerosol contamination of the data. 

C. THE MODEL 

LOWTRAN7 (Low Resolution Transmittance) is a model that calculates 
atmospheric transmittance, scattering and absorption. It incorporates separate 
molecular profiles for all major as well as 13 minor and trace gases, all as a 
function of altitude (0 to 100km). Representative atmospheric aerosol, cloud, 
and rain models are provided with the ability to replace any or all values with 
theoretical or measured data at the discretion of the user. The model also 
provides six reference atmospheres with options to redefine or retain profiles. 

LOWTRAN7 is a band model. Average transmittance is calculated over 20 
cm'l bands. The band parameters have been determined semi-empirically 



3 



determined semi-empirically through the averaging of line parameters. In 
principle, the line by line method is available, but the calculations are too 
elaborate and time consuming for most users. Some accuracy is gained by a line 
method, but it is not significant enough to warrant the loss of efficiency. 
Individual lines within a band may vary significantly, yet it is fair to characterize 
the effective band parameters through the averaging of line parameters so long 
as a specific line result is not desired. 

For use in multiwavelength applications (ie. satellites) the appropriate 
spectral response functions for the sensor can be applied to the output data 
through LOWFIL7 (LOWTRAN7's filter program). The data set derived 
through the running of the main program is filtered through the sensor's 
response function and produces a normalized radiance estimate of the sensor 
measurement. 
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II. THEORY 



A. PARTICLE SIZE INDEX (S 12 ) 

The Advanced Very High Resolution Radiometer (AVHRR) has been a part 
of all the NOAA polar orbiting satellites. It has five channels: 



TABLE 1 

NOAA-7 AVHRR CHANNEL BANDWIDTHS 



Channel # 
1 

2 

3 

4 

5 



Bandwidth (urn) 
0.58 - 0.68 
0.725 - 1.10 
3.55 - 3.93 
10.30- 11.30 
11.50- 12.50 



AVHRR satellite detected radiances in the red-visible (chan. 1) and the near 
infra-red (chan. 2) wavelength ranges have been found to be positively 
correlated with aerosol optical depth (8) (Durkee 1984; Griggs 1979). 

Using a single scattering approximation and knowing that optical depths over 
the ocean are small (Pfeil 1986), Frost (1988) used a simplified form of the 
Radiative Transfer Equation (RTE): 

MS,h, 4>) ' 5 ^jp L P<©) 8 (2.1) 

where, L = diffuse radiance, 

(Do = single scatter albedo. 
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F 0 = solar radiative flux, 

fi = cos0 (0 = satellite zenith angle), 

p = scattering phase function, 

© = single scattering angle, 

8 = optical depth. 

In an optically thin (small 8) atmosphere, multiple scattering effects are very 
small, making single scattering a reasonable assumption (Frost 1988). 

There are three primary sources of upwelled radiance (L) that contribute to 
satellite detected radiance in the red-visible and near-infrared wavelengths. The 
contributions can be expressed as: 

L = La + Lp + Ls (2.2) 

where L a is aerosol scatter, L r is molecular or Rayleigh scatter, the wavelength 
of radiation is much greater than the size of the scattering material, and Ls is sea 
surface reflectance. Except for sunglint, Ls is very small for red wavelengths 
(albedo = 0.5%) and zero for wavelengths greater than 0.7|im (Ramsey 1968). 
Therefore, as long as sunglint geometry is avoided, total radiance is due to the 
sum of aerosol and molecular scattering effects (L = L a + Lp). 

L r is strongly wavelength (X) dependant, but does not vary spatially. 
Therefore L r forms a baseline L(X) and aerosol particles are the dominant 
source of upwelled radiance variations (AL ~ AL a ). From Eqn. 2.1, changes 
in satellite measured radiance are directly determined from changes in aerosol 
characteristics like too, p(0), and aerosol optical depth (8 a ). Therefore, to 
retrieve a more accurate 8 from L a measurements, knowledge of p(0) 
variations is required. 
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Durkee (1984) showed that retrieved radiances at red wavelengths (chan. 1) 
are more sensitive to varying densities of small aerosol particles than that of 
near-infrared (nir) wavelengths (chan. 2). The most significant parameter 
contributing to retrieved radiances at AVHRR channel 1 and 2 wavelengths is 
the extinction coefficient (a e xt) and specifically the vertical integral of a e xt or 
the optical depth (8). 



(Text is the sum of extinction due to absorption (a a bs) and due to scattering 
(reseat)- For marine aerosols, absorption is negligible, therefore aext ~ ctscat 
and 



H 




(2.3) 



H 




(2.4) 



Three terms contribute to Gscaf 



CO 




(2.5) 



where, tct 2 = particle cross sectional area, 



Qscat(m,X,r) = scattering efficiency, 
m = complex index of refraction, 
X = wavelength, 
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r = particle radius. 



n(r) = number distribution of particles. 

Figure 1 shows the influence of the terms of Eqn. 2.5 as a function of r and 
the resulting number-weighted extinction (rcr 2 -Q S cat‘ n ( r )) f° r channel 1 and 
channel 2. Scattering cross section (tir 2 -Qscat) and the particle size distribution 
(n(r)) for a modeled distribution of marine particles at 80% RH are plotted. As 
wavelength increases, so does the r of most efficient scattering (Qscat)- The 
result is a shift in the scattering cross section and a similar shift in the r that most 
significantly contributes to a S cat, and therefore L at a given X. Therefore, the 
slope of L(X) closely resembles the slope of S(X), and is similar to the slope of 
the aerosol particle size distribution (n(r)). Since aerosol particles most 
effectively scatter radiation at a wavelength similar to the radius of the particle, 
an increase in the number of small aerosol particles will increase the ratio of 
L c hl to L c h2> and the spectral variation can be quantified: 

tflo'Fp - 

L(red) _ L(chl) _ 4n P< e > 8 (chQ _ 5^ (2 fi) 

L(nir) L (ch2 ) ^ p(e)5(ch2) »G*> 



This ratio is called the Particle Size Index (S 12 ) (Frost 1988): 1 

<5 _ (L)red _ (5)red 

^ (L)nir (5)nir 



(2.7) 



tWhile Frost removed Lr prior to taking the ratio, it will not be removed in 
this study. It shall be accounted for in the parameterization of p(@). 
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Fig. 1 Extinction formulation as a function of particle radius for 
channels 1 (solid) and 2 (dash). Modified from Durkee (1984). 
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f _UJD) (J) 



Satellite measured radiances therefore cannot provide the actual particle size 
distribution, but the S 12 gives the slope of the size distribution curve. In 
general, the greater the contribution of smaller particles, the larger the ratio. 

Figure 2 shows a plot of NOAA-7 AVHRR channel 1 and channel 2 
response function curves superimposed on LOWTRAN7’s Tropical model 
atmosphere output radiances (L(X)). Although as this paper will show, water 
vapor absorption does effect S 12 , the dominant variable influencing L(^.) in 
these wavelength bands is aerosol scattering. For the red and near-infrared 
region of the spectrum, variances in L(X) are proportional to changes in 6 a (X) 
and therefore aerosol scattering. Note that L(X) in the channel 1 window (short 
wavelengths, small particle scattering) is greater than that of channel 2. This is 
primarily due to the slope of L r (X) and leads to an S 12 ratio that is generally 
greater than 1 . 

B. WATER VAPOR ABSORPTION IN THE VISIBLE / NEAR 

INFRA-RED 

Although aerosol scattering is the dominant influence on satellite retrieved 
radiances at visible and near-infrared wavelengths, the role of atmospheric water 
vapor is not insignificant. The magnitudes of the overall L in these wavelengths 
are not appreciably effected, but often the goal of using satellites to estimate 
aerosol particle characteristics is to look for variations, which are an order of 
magnitude smaller than the total L. 

In order to observe sometimes subtle, yet often significant trends such as 
global climate variations, an improvement in the accuracy of aerosol estimations 
is necessary. 
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Fig. 2 NOAA-7 AVHRR channel 1 and 2 response functions and total 

atmospheric radiance. 

Figure 3 shows the effect of water vapor alone on x for AVHRR channels 1 
and 2. Notice that for channel 1, x is greater than channel 2, meaning that less 
absorption is taking place. This would tend to decrease L2 more than Li and 
thus lead to an over-estimate of the S12. Therefore a change in water vapor in a 
region may be misconstrued as a change in aerosol characteristics or disguise a 
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significant change. This implies that we might expect to see a minimum S[ 2 (w) 
at w = zero, and an increase in Si 2 (w) with w. 




Wavelength (microns) 



Fig. 3 NOAA-7 AVHRR channel 1 and 2 response functions and 
atmospheric transmittance due to water vapor. 

In Figure 4, the influence of water vapor is again demonstrated when water 
vapor is added to LOWTRAN7’s Tropical atmosphere model. Notice that 
while radiances decrease in both channel wavebands, channel 2 is more 
significantly effected. This again will cause Sn to increase with w even though 
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Li and L 2 decrease. There does appear to be some increase in L below 0.6 pm. 
This is most likely due to an effective shift in particle size due to coalescence of 
water vapor. This study intends to look at the influence of water vapor alone on 
retrieved radiances, and aerosol growth will be inhibited. 
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Fig. 4 NOAA-7 AVHRR channel 1 and 2 response functions and the 
influence of water vapor on atmospheric radiance. 



C. BRIGHTNESS TEMPERATURE DIFFERENCE (T 4 -T 5 ) 

In order to take advantage of the water vapor parameterization of the 
scattering phase function proposed in this study, some method of measuring the 
water vapor present would be required. 

A significant benefit of using the AVHRR is that while channels 1 and 2 
provide information on optical depth variations, channels 4 and 5 in the 10 to 13 
pm window region have demonstrated usefulness in atmospheric water vapor 
retrieval. 

Unlike the aerosol backscattering of channels 1 and 2, the brightness 
temperature is the result of single direction radiation transfer, from the sea 
surface to the sensor. Therefore, T is directly proportional to retrieved T, with 
absorption reducing the radiance reaching the sensor. 

From Dalu, et al. (1981): 

T s = Tcc + g(w)(T a -Tp) (2.8) 

T s = Sea Surface Temperature (SST) 

T<x = Brightness Temp, of most transparent channel 
Tp = Brightness Temp, of second channel 
g(w) = function of water vapor content 
1 - Ta 

= [ c (1 - Tp) - (1 - Ta) 1 

where Ta and Tfl are total atmospheric transmittances for channels a and (3 and 
C=(T S - T~p,)/(Ts- fa)- Solving for(T a -Tp): 

rr Tn\ ( T s- T a) r ~ 

(T a -Tp)= g(w) (2.9) 
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By using an approximation proposed by Prabhakara, et <a/.( 1 974): 

(Ts-Tcc) = (T s - T a )Kw sec0 (2. 10) 

where 0 is the scanning angle and K is an absorption coefficient. Column water 
vapor (w = Jp w dz) (g m'2) can be related to the brightness temperature 
difference: 



(Ts- T a) 

(T a -Tp)=K g(w) w sec0 (2.11) 

Although both g(w) and (T s - T a ) increase with increased water vapor, their 
ratio remains nearly constant (Dalu 1986). This leaves us with: 

w = A(Ta-Tp)cos0 (2.12) 

where A = const. = (g(w)/(K(T s - T «))• 

To maximize the effectiveness of this scheme, T a should be minimally 

influenced by water vapor in the atmosphere to allow the retrieval of water 
vapor temperatures near the surface. Tp, on the other hand, should be more 

sensitive to the atmospheric water vapor distribution. Under such conditions, 
(T<x-Tp) would indicate the influence of water vapor, which would be 

proportional to the amount (w). 

Channel 4 (~ll|j.m) and channel 5 (~12|im) windows of the AVHRR were 
designed for sensing of the sea surface temperature (SST) corrected for 
atmospheric absorption (Dalu 1986). Despite being optimized for SST retrieval, 
they describe a linear function from which the atmospheric water vapor 
correction can be derived (McMillin and Crosby 1984). Figure 5 shows 
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NOAA-7 AVHRR channels 4 and 5 response function curves superimposed on 
three model atmosphere x of varying water vapor amounts. Notice that T4 is 
generally greater than 15, leading to a positive T 4 -T 5 , and that higher T 4 -T 5 



correlates directly with higher water vapor densities (Tropical > US Standard > 
Midlatitude Winter). 




Fig. 5 NOAA-7 AVHRR channel 4 and 5 response functions and 
atmospheric transmittance due to water vapor. 
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Although not the ideal sensor windows for water vapor retrieval (the 
channels are too similar), the use of T4-T5 is justified to within the accuracy of 

the approximations made to achieve Eqn. 2.12. The most beneficial aspect of 
T4-T5 is the co-location and temporal linkage with channels 1 and 2, providing a 

water vapor snapshot simultaneously and using the same path as the aerosol 
optical depth (5) measurements. 

Therefore, an estimate of vertically integrated water vapor density (w) is 
provide by: 



w = A*(T4-T5 )cos0 (2.13) 

Using A* = 19600 g /(°K m^) and a scanning angle 0 = 0, Dalu (1986) applied 
this relationship to a radiative transfer model using a variety of temperature and 
humidity profiles. Figure 6 shows the correlation coefficient for the resulting 
line is R = 0.99 with the error given as ± 1.5 kg m'2. Comparing results 
calculated using radiosonde data with shipboard measured water vapor amounts, 
Figure 7 shows the correlation coefficient dropped to R = 0.78 with an error of 
± 4.0 kg m‘2. The decrease in accuracy most likely resulting from calibration 
errors as well as inexact temporal and spatial correlation of ship and satellite 
data. 

At channel 4 and 5 wavelengths, scattering is a relatively insignificant 
influence determining brightness temperature (T). The dominant factor is water 
vapor absorption of surface temperatures. 
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Fig. 6 Total water vapor content (kg/m 2 ) as a function of brightness 
temperature difference from Dalu (1986). 

Some of the long wave radiation emitted by the sea surface is absorbed by 
water vapor. The vapor then re-emits radiation based on temperature in 
accordance with the Planck function. If the atmosphere, and the water vapor 
included, were not generally cooler than the sea surface, the absorbed radiation 
would be re-emitted at the same rate, not having any effect on radiance or 
brightness temperature. 
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Fig. 7 Comparison of water vapor content from ships with radiative 
transfer model runs using ship radiosonde data. From Dalu (1986). 
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III. PROCEDURE 



A. INTRODUCTION 

This study uses the LOWTRAN7 model to improve algorithms used in 
satellite retrieved optical depth measurements, specifically a correction due to 
varying amounts of water vapor present. The benefit in using a theoretical 
model is that the parameters involved can be specifically quantified and 
controlled. Therefore, a precise relationship between satellite measured 
radiances (L, w/cm 2 -steradian-cm' 1 ) and varied parameters such as optical depth 
(8), atmospheric temperature (T, °K), column water vapor (w, gm m' 2 ) can be 
determined. These relationships will be consistent with that of the real 
atmosphere to the degree that the model represents the real atmosphere. 

B. OVERVIEW 

Figure 8 is a flow chart of the procedure beginning with a user defined 
atmosphere of known quantities of aerosol optical depth (5), water vapor 
densities (p w ), and surface temperatures (T S urf)- Three programs are run for 
each data set. Both LOWTRAN7 and it's Filter function program, LOWFIL7, 
make use of formatted Files (TAPEY.DAT) for input, output, and to pass coded 
information between the programs. The program RETRIEVE was specifically 
designed for this study, to compile data in a condensed format. 

LOWTRAN7 is supplied a modeled atmosphere, sun/satellite/earth 
geometry, as well as a desired wavenumber (k) band. Wavenumber dependent 
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TAPE6.DAT 

Broad Band Radiances 



TAPE5.DAT 

► 850 km satellite height 

► 10 km troposphere 

► Background stratospheric extinction 

► User defined aerosol scattering phase function 

► Varied column water vapor (ill) 

* Varied , user defined aerosol extinction 
coefficients (xlO) 

► 260°K tropospheric temperature 

► Varied surface temperature (x4) 

► NOAA-7 AVHRR frequency windows (x2) 

Chan. 1.2 
Chan. 3.4 



I 



LOWTRAN7 



TAPE5.DAT 

NOAA-7 AVHRR 
Chan. 1,2 
Chan. 4,5 

Sensor Response Functions 



TAPE7.DAT 

Broad Band Filter Input 






LOWFIL7A/B 



T 



TAPE6.DAT 

Filtered Radiances: 

1 )Total Integrated 
2)FUtcr Normalized 




RETRIEVE.DAT 

Chan. 1 Radiance = 

Chan. 2 Radiance = 

L1/L2 = S 1 2 



LI 
L 2 



Chan. 4 Brightness Temp = T4 

Chan. 5 Brightness Temp = T 5 

T4-T5 = T4-T5 



Fig. 8 Flow chart of model run. 
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radiances are calculated considering scattered, absorbed, and emitted radiation as 
it would be "seen" by a perfect sensor at the satellite location. Figure 9 shows 
the AVHRR sensor response function curves (<D(k)) that LOWFIL7 combines 
with coded information derived from the previous LOWTRAN7 run to 
calculate the effective filtered radiances as a function of wavenumber or 
wavelength. The sensor weighted total integrated radiance is then determined. 

Ltot = |L(K)<t>(K)dK (3.1) 

LOWFIL7 was modified for use in this study to calculate the integral of the 
sensor's response function which is then used to calculate a normalized radiance 
(LOWFIL7A/B). 2 



jL(K)<l>(K)dK 

L(chan) - f (3.2) 

J<P(K)dK 

The program RETRIEVE reads the output files from LOWFIL7A/B and 
writes only the normalized radiances for channels 1, 2, 4, and 5 (Li,L 2,L4,L5). 
The ratio of channel 1 to channel 2 radiances is calculated to determine the 
particle size index parameter (S12) (Frost 1988). 3 



q _ (L)red _ El 

bU " (LW ' E 2 



(3.3) 



2 LOWFIL7 normalizes by , <Pmax instead of J<b(K)dK. 

3 Unlike Frost, this study does not first remove L r . 
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Relative Response Relative 




Wavelength (microns) 




Fig. 9 NOAA-7 AVHRR filter function response curves. 
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Channel 4 and channel 5 radiances are converted to brightness temperatures 

(T b ). 



T b(4,5) = 3 ( 3 - 4 ) 

In ( 1 + ) 

L(chan 4,5) 

Where Tb(4,5) is the brightness Temperature (K) for the radiance L(4 t 5) 
(w/cm^-steradian-cm* 1 ), k is the central wavenumber of the channel filter 
( K chan4 = 927.22 cm* 1 and Kchan5 = 840.872 cm* 1 ), and Cl and C2 are constants 
(Ci= 1.1910659x10* 1 2 w/cm^-steradian-cm* 4 and C2= 1.438833 cm* 1 K) 
(NOAA Polar Orbiter Data 1986). Finally, the brightness temperature 
difference between Tb4 and Tb5 is calculated (T4-T5). 

C. INPUT CHARACTERISTICS 

The modeled atmospheres were restricted to cloud-free marine environments 
where backscattered radiance is dominated by aerosol particles (Durkee 1984). 
Since the satellite senses the effect of vertically integrated parameters, the models 
used were homogeneous from the surface to the tropopause. For ease of 
integration, the model tropopause was set to 10 kilometers (the average height of 
the mid-latitude tropopause is about 11 kilometers (Stull 1988)). From the 
tropopause to the satellite height, only background extinction due to meteoric 
dust is represented from 30 to 100 kilometers. The following are the important 



user-defined parameters: 

aerosol optical depth (5a(A.)) varied 

column water vapor (w) varied 

tropospheric air temperature (Tair) 260 K 

surface temperature (T sur f) varied 
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surface albedo (Salb) 0.0 (blackbody) 

satellite height (H) 850 km 

season spring-summer 

marine phase function (P(©)) from Shettle and Fenn (1979) 

wavenumber windows, resolution 

chan 1,2 8400-20050 cm* 1 , 50 cm* 1 

chan 4,5 770-1005 cnrl, 5 cirri 

Molecular constituents such as O3, CH4, N2O, CO, CO2, O2, NO, SO2, NO2, 
NH3, and HNO3, used LOWTRAN7 model default profiles for midlatitude 
summer and were kept constant for all model runs. The sun/satellite geometry 
remained the same throughout all the model runs with the satellite and sun in a 



straight line perpendicular to the earth's surface. 

1. Optical Depth Profiles 

Satellite windows in the red-visible and near-infrared sense the result of 



vertically integrated parameters. The most significant influence on retrieved 
radiances is the vertical integral of (Text or $. To study this effect, ten profiles 
of 8 with constant a e xt to 10 km were designed. 

Using the relationship 



CTscat ** (3.5) 

and an Gscat of .05 km -1 at A, = .55 |im, a profile of Gscat as a function of X was 
derived (a S cat(X)). After integrating over the 10 km height to determine 8, 
Figure 10 shows this curve (M2) and four multiples of this profile that complete 
the first set of optical depth profiles. 
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MULTIPLE PROFILE 




Fig. 10 Multiple optical depth profiles (M). 

Figure 1 1 shows a second set of profiles created using the same relationship and 
adding equal quantities to create parallel profiles. 

2. Water Vapor 

All the water vapor was contained in the homogeneous troposphere (0 
to 10 km). The maximum value was selected from the tropical atmosphere 
model included in LOWTRAN7 (3.322 gm cm*2). User-defined profiles are 
input as water vapor mass density (pw) (gm nr 3), which after integrating in the 
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PARALLEL PROFfLE 




Fig. 11 Parallel optical depth profiles (P). 
vertical gives column water vapor (w) of equivalent magnitudes in units of gm 
cm'2. 

H 

w = Jp w dz (3.6) 

Ten intermediate quantities between 0 and 3.322 gm cm*2 were selected with 0 
included as the eleventh. 
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3. Surface Temperature 

T4-T5 is a measure of water vapor (Jacobowitz and Gruber 1990) and 
will be compared with the user input water vapor concentrations, as an empirical 
method of estimating w. T5 is more sensitive to water vapor effects on 
atmospheric transmittance while T4 senses more of the near-blackbody, long- 
wave radiation from the surface. In order to demonstrate relationships, the 
model required a warmer surface temperature than the atmosphere when the 
troposphere is isothermal. With equal atmospheric temperature (T a tm)and 
surface temperature (T sur f), the measured brightness temperatures would be 
identical. Also, surface temperatures colder than the temperature at 10 km are 
unrealistic. Therefore, while the tropospheric temperature was held constant, 
the surface temperature was varied to achieve four representative modeled 
Tatm/Tsurf relationships: 

Tatm 260°K 

T sur f(w) 270°K 

T sur f(x) 280°K 

T SU rf(y) 290°K 

T sur f(z) 300°K 

D. PROCEDURAL SUMMARY 

The model was run at the Naval Postgraduate School's (NPS) Interactive 
Environmental Digital Analysis Laboratory (IDEA LAB). A complete set of 
input files consisted of two sets of five 8 profiles, with eleven w variations, four 
surface temperature conditions, and two wavenumber windows. These 880 files 
(2X5X 11 X4X2) generated 1760 intermediate and final output files, with 
the significant data 'retrieved' and collated in one resultant array. All the 
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calculations were done as a large batch job requiring about 4.5 hours to run. 
The job was completed in both the single and multiple aerosol scattering modes 
of LOWTRAN7. 
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IV. RESULTS 



A. THE EFFECT OF AEROSOL OPTICAL DEPTH VARIATIONS 
ON S12 

In general, an increase in wavelength dependant aerosol optical depth 
(§a(^)) resulted in an increase in Li and L2 while S12 values decreased. Since 
S12 is generally greater than one (Li>L 2). Increasing the radiance in both 
sensor windows by similar magnitudes has a greater influence on the smaller L, 
driving S12 toward one. This relationship will remain true because of the 
Rayleigh contribution, which is greater for Lj. 

1. Multiple 6 a (^) Scheme 

Increasing 5a profiles result in decreasing S12 ratios, with the rate of 
decrease diminishing due to the effect of increasing 5 a (X) slopes. Figure 12 
shows the S12 variations resulting from the multiple 5a profiles of Fig. 10. 

2. Parallel 5 a (X) Scheme 

By constraining the curves to be parallel, an increase in 8a is an equal 
increase at all wavelengths. The S12 ratios decrease at a more rapid rate than in 
the previous scheme since a change in slope is no longer limiting the decrease in 
S12. Figure 13 shows the S12 variations resulting from the parallel 5a profiles 
of Fig. 11. 

B. THE EFFECT OF SURFACE TEMPERATURE VARIATIONS ON 

THE T4-T5 REPRESENTATION OF W 

Although the modeling accomplished in this study was not designed to prove 
Dalu's technique for determining vertically integrated water vapor density (w) 
from T4-T5, some significant aspects of his theory are supported by the results. 
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Fig. 12 S12 versus column water vapor (w) for the multiple optical 

depth scheme (M). 

Figure 14 shows the T4-T5 values derived from modeled atmospheres of known 
w given a 260K isothermal air temperature (T a ir) and various T$urf conditions. 
The slopes of the T4-T5 varied from Dalu's estimate, yet the linear relationship 
of Eqn. 2.13 seemed evident. 

The variation in the slopes is a result of the assumptions made in the model 
atmospheres. The use of an isothermal 10 km atmosphere is reasonable for 
NOAA-7 AVHRR channel 1 and 2 retrievals since these wavelengths 
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Fig. 13 S 12 versus column water vapor (w) for the parallel optical 

depth scheme (P). 

(red-vis, nir) sense the vertically integrated aerosol backscattering with little 
height dependance. Since absorption plays a big role in channel 4 and 5 
brightness temperature retrievals in the actual atmosphere, Tb 4 and Tb 5 results 
are significantly influenced by vertical temperature variations. The isothermal 
assumption used in this model study, therefore, causes quantitative T 4 -T 5 results 
to be unrepresentative of the real atmosphere. However, qualitative results such 
as the linear relationship are retained. Table 2 compares the slopes derived from 
applying Dalu’s theory to each T a ii7Tsurf variation. 
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Fig. 14 T 4 -T 5 versus column water vapor (w) for a variety of 

surface temperatures. 
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TABLE 2 

MODEL OUTPUT COEFFICIENTS APPLIED TO DALU'S 
THEORY (w = A(T 4 -T5)) 



Tair/Tsurf (K/K) 



A (g K-l m-1) 



260/300 



.7615 



260/290 



.9944 



260/280 



1.4995 



260/270 



3.6653 



Dalu 



1.9600 



C. THE W CORRECTION TO S12 

In order to present the correction to satellite retrieved S12 values required 
due to the presence of water vapor, the results will be displayed first graphically, 
then in a new parameterization that can be applied to the aerosol scattering phase 
functions used in retrieval of aerosol characteristics from satellites. 

1. Graphical Application 

RETRIEVE lists the following results for each modeled atmosphere: 



A set of output variables was calculated for each variation of input parameters. 
The combined results are displayed in Fig. 15. 



Input Parameters: 



Output Variables: 



5 a (X) 

Tsurf 




w 
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Both the multiple and parallel S 12 results are displayed and seem to follow a 
similar relationship to w. Note that M2 and PI curves are coincident. They use 
the same input 8 a (A.) values that formed the basis for each set of curves. 

2 . The Reformulated Parameterization 

Looking at the S 12 values plotted for varying w and 8 a (X) curves, 
there appears to be a functional relationship. It is apparent that the water vapor 
influence of S 12 is not dependent on the shape of the optical depth profile, but 
was a function of w and the magnitude of S12 when w = 0 (Si 2 (dry)). Figure 
16 shows that when the ratio of Si 2 :Si 2 (dry) was plotted as a function of w, the 
shape of the curve was conserved for all model outputs (multiple and parallel). 
The slight variation at higher w is partly due to computational round off errors. 
The curve that best fits the data is: 

Sn(dry) = (1 + .0332 Vw) S 12 (4.1) 

with r^ = 0.995. 

From Eqn. 4.1, measured S 12 and estimated w from T 4 -T 5 determines 
Si 2 (dry). 

D. SINGLE VERSUS MULTIPLE SCATTERING 

All data sets were run with both single scatter and multiple scatter 
approximations. The results were calculated to the seventh decimal place and 
were precisely the same for either case. The model, therefore, supports Frost's 
use of a simplified form of the Radiative Transfer Equation (RTE) that assumes 
single scatter (Eqn. 2.1). 
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*2 

Column Water Vapor (g cm ) 



Fig. 16 S12 normalized by Si2(dry). 
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V. CONCLUSIONS 



Present techniques for satellite retrieval of aerosol characteristics using the 
AVHRR carried on the NO A A series platforms require that analyzed pixels be 
free of clouds. This is reasonable since liquid water effects on radiative 
processes are quite variable and difficult to measure. Clouds at any level in a 
satellite scene can obscure the sources it was designed to sense, or saturate the 
sensor with reflected radiance. It is understood that water in a gaseous form 
would still be present, but its influence on L in the wavebands most effectively 
used to estimate aerosol parameters (Chan. 1 and Chan. 2) has been assumed to 
be insignificant in previous studies. 

Although aerosol scattering is still the dominant influence on L in these 
wavebands, this study has shown that water vapor does introduce an error, and 
that the error can be removed. This does not invalidate the present techniques, 
but provides a significant improvement that can be applied through the use of the 
variable scattering phase function introduced by Frost (1988). Since aerosol 
variations can be of global scale, even very small changes may have global 
ramifications. 

Through the use of the LOWTRAN7 atmospheric propagation model, this 
study has shown support of the following techniques and theories presently 
accepted: 

• The concepts of Frost's particle size index parameter (S 12 ) using 
AVHRR channels 1 and 2. 

• Frost’s use of a simplified radiative transfer equation based upon an 
assumption of single scatter. 

• Dalu's water vapor retrieval technique using AVHRR channels 4 and 5 
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By introducing varying quantities of water vapor to the modeled 
atmospheres of differing optical depth (5) and T sur f/T a j r profiles, additional 
information has been generated: 

• The error introduced to S 1 2 due to water vapor is demonstrated and 
quantified. 

• The error has been shown to be a function of Si2(dry) and column 
water vapor (w). 

The influence of optical depth variations can be isolated from the variation 
due to w through the use of the variable scattering phase function (P(©)) 
proposed by Frost, (1988). Using Eqn. 2.1, S 12 with 2.0gm cm' 2 of water 
vapor produces a 15% error in P(@) and therefore retrieved aerosol optical 
depth (8 a ). 

By using a model to determine water vapor's influence on satellite-measured 
aerosol characteristics, the parameters normally estimated can be strictly 
controlled. As known quantities, their theoretical relationships to satellite- 
measured quantities can be observed. The next step is to apply the correction 
discovered here to empirical satellite data with verifying in situ data sources such 
as aircraft or ships. 



39 



LIST OF REFERENCES 



Bloembergen, N., Patel, C. K., and G. Pake (Chairpersons), 1987: Science and 
Technology of Directed Energy Weapons. Report of the American 
Physical Society Study Group, Reviews of Modern Physics, 59, S1-S202. 

Charlson, R.J., Lovelock, J.E., Andrea, M.O., and S.G. Warren, 1987: Oceanic 
Phytoplankton, atmospheric sulfur, cloud albedo and climate. Nature, 
325, 655-661. 

Dalu, G., 1986: Satellite Remote Sensing of Atmospheric Water Vapor. 

International Journal of Remote Sensing, 7, 1089-1097. 

Dalu, G., Prabhakara, C., and Lo, R. C., 1981: Improved accuracy of the 
remote sensing of sea surface temperature. Oceanography from Space , J. 
F. R. Gower, ed.. New York:Plenum, 109. 

Durkee, P.A., 1984: The relationship between marine aerosol particles and 
satellite -detected radiance. Ph.D. dissertation, Colorado State University, 
Fort Collins, CO, US ISSN 0067-0340, 124pp. 

Frost, E.M., 1988: Global scale estimates of aerosol particle characteristics. 
M.S. Thesis, Naval Postgraduate School, Monterey, CA., 54pp. 

Griggs, M., 1979: Satellite Observations of Atmospheric Aerosols During the 
EOMET Cruise. Journal of Atmospheric Science., 36, 695-698. 

Jacobowitz, H., and A. Gruber, 1990: Estimation of the Outgoing Longwave 
Flus from NOAA AVHRR Satellite Observations. Seventh Conference on 
Atmospheric Radiation. American Meteorological Society, 4pp. 

McMillin, L. M., and D. S. Crosby, 1984: Theory and Validation of the 
Multiple Window Sea Surface Temperature Technique. Journal of 
Geophysical Research, 89, 3655-3661. 

NOAA POLAR ORBITER DATA (TIROS-N, NOAA-7, NOAA-8, NOAA-9 
AND NOAA-10) Users Guide. U. S. Department of Commerce, National 
Oceanic and Atmospheric Administration, December 1986. 



40 



Pfeil, F. R., 1 986: Developing a physical basis for an aerosol climatology of the 
Pacific Ocean. M.S. Thesis, Naval Postgraduate School, Monterey, CA, 
76pp. 

Prabhakara, C., Dalu, G., and Kunde, V. G., 1974: Estimation of sea surface 
temperature from remote sensing in the 11 to 13 pm window region. 
Journal of Geophysical Research, 79,5039. 

Ramsey, R. C., 1968: Study of the Remote Measurement of Ocean Color. Final 
Report, TRW, NASW-1658. 

Shettle, E. P., and R. W. Fenn, 1979: Models for the Aerosols of the lower 
Atmosphere and the Effects of Humidity Variations on Their Optical 
Properties. AFG-TR-79-0214 Air Force Geophysics Laboratories, 
Hanscom AFB, MA., 94pp. 

Stull, R. B., 1988: An Introduction to Boundary Layer Meteorology. Kluwer 
Academic Publishers, Dordrecht, The Netherlands, 666 pp. 



41 



INITIAL DISTRIBUTION LIST 



No. Copies 

1 . Defense Technical Information Center 2 

Cameron Station 

Alexandria, V A 22314 

2. Library, Code 52 2 

Naval Postgraduate School 

Monterey, CA 93943-5000 

3. Chairman (Code OC/Co) 1 

Department of Oceanography 

Naval Postgraduate School 
Monterey, CA 93943-5000 

4. Chairman (Code MR/Hy) 1 

Department of Meteorology 

Naval Postgraduate School 
Monterey, CA 93943-5000 

5. Professor P. A. Durkee (Code MR/De) 1 

Department of Meteorology 

Naval Postgraduate School 
Monterey, CA 93943-5000 

6. Professor K. L. Davidson (Code MR/DS) 1 

Department of Meteorology 

Naval Postgraduate School 
Monterey, CA 93943-5000 



7. LT Timothy P. Mahony, USN 1 

NAVOCEANCOMCEN Box 31 
FPO New York, NY 09549-3200 



42 



1 



8. Commanding Officer 

Naval Oceanographic and Atmospheric 
Research Laboratory 
Stennis Space Center 
MS 39529-5004 



9. Director 1 

Naval Oceanographic and Atmospheric 
Research Laboratory 
Monterey, CA 93943-5006 



10. Chairman 1 

Oceanography Department 
U. S. Naval Academy 
Annapolis, MD 21402 



1 1 . Office of Naval Research (Code 420) 1 

Naval Ocean Research and Development 
Activity 

800 N. Quincy Street 
Arlington, VA 22217 



12. Scientific Liaison Office 1 

Office of Naval Research 
Scripps Institution of Oceanography 
La Jolla, CA 92037 



13. NOAA Library 1 

7600 Sand Point Way NE 
Building 3 
Seattle, WA 98115 



43 




) /(cTr 




r -e-b I ( I ^ ~i 



Thesis 

M277175 Mahoney 

o ^ 1 TJ*nt-o>- y o n n »■ '' 

on satellite-measured 
aerosol characteristics . 



